Phenomenological study of the electronic transport coefficients of graphene 
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Using a semi-classical approach and input from experiments on the conductivity of graphene, we 
determine the electronic density dependence of the electronic transport coefficients - conductivity, 
thermal conductivity and thermopower - of doped graphene. Also the electronic density dependence 
of the optical conductivity is obtained. Finally we show that the classical Hall effect (low field) in 
graphene has the same form as for the independent electron case, characterized by a parabolic 
dispersion, as long as the relaxation time is proportional to the momentum. 
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Introduction. Since the experimental measurement of 
an electric field effect in graphene^ and the observation 
of the odd integer quantum Hall effect^i 3 - the electronic 
properties of graphene have been attracting a consider- 
able attention from the community. Recent qualitative 
reviews^^ on the physics of graphene give a brief ac- 
count on both the experimental and theoretical status of 
the field. 

The measured conductivity of graphene has two dis- 
tinct fingerprints: a linear dependence upon the gate 
voltage of the DC conductivity, and a minimum con- 
ductivity at the neutrality point of one electron per car- 
bon atom. The value of the conductivity minimum is 
of the order of the quantum of conductance e 2 /h. This 
minimum of conductivity is a quantum mechanical ef- 
fect which comes about as a consequence of disorder 7 . 
Disorder promotes a finite density of states at the Dirac 
point which is responsible for the minimum of conduc- 
tivity, a mm- Unfortunately there are a number of dif- 
ferent theoretical predictions available for cr m in in the 
literature^^^ 2 -! 3 -^^^. Except for the numerical 
work of Nomura and MacDonaldi 7 - all the calculations 
for the conductivity minimum are off by some numerical 
factor from the experimentally measured value. 

Also interesting by itself is the linear dependence of 
the conductivity on the gate potential. Since the gate 
potential depends linearly on the electronic density, n, 
one has a conductivity er oc n. As shown by Shon and 
Andoi£ if the scatterers are short range one obtains a 
DC conductivity that is independent of the electronic 
density, at odds with the experimental result. Nomura 
and MacDonal d 17 ' 19 showed that considering a scatter- 
ing mechanics based on screened charged impurities it is 
possible to obtain from a Boltzmann equation approach 
a conductivity varying linearly with the density, in agree- 
ment with the experimental result. In this brief report 
we address the problem of finding the electronic density 
dependence of several transport properties of graphene. 
Since the electronic density is easily controlled by a gate 
voltage the expressions we derive here can certainly be 
tested experimentally. We therefore present in what fol- 
lows a calculation of the electronic density dependences 



of the DC conductivity, the thermal conductivity and the 
thermopower. In addition, we study how the electronic 
density enters in the optical conductivity. Finally, we 
show that a certain form of the scattering rate is neces- 
sary if the classical expression for the Hall effect is to be 
maintained, in agreement with the experiments. A sim- 
ilar semi-classical approach was used by Falkovsky 2 ^ to 
study the temperature dependence of the Hall conduc- 
tivity in graphene. The link between the Kubo-Streda 
formulation and the semi-classical Boltzmann equation 
approach was discussed, in the context of massive Dirac 
fermions, by Sinitsyn et al.Q 

Boltzmann equation and relaxation times. The Boltz- 
mann equation has the form 22 

-v k -V r .f(e k )-e/h(E+v k xHyV k f(e k ) = - ^ 

scatt. 

_ (1) 

The solution of the Boltzmann equation in its general 
form is difficult and one needs therefore to rely upon 
some approximation. The first step in the usual approx- 
imation scheme is to write the distribution as f(e k ) — 
/°(cfc) + ff(efc) where /°(efc) is the steady state distribu- 
tion function and g(e k ) is assumed to be small. Inserting 
this ansatz in Eq. ([1]) and keeping only terms that are 
linear in the external fields one obtains the linearized 
Boltzmann equation 22 which reads 



df°(e k ) 
de k 

dfh 
dt 



v k 



e -^\ V r T + e( E -lv r C 



v k ■ V r .g fc + ~(v k X H) ■ Vfcgfc . (2) 

n 



The second approximation has to do with the form of the 
scattering term. The simplest approach is to introduce 
a relaxation time into the formalism. This is done by 
considering the approximation 



dfh 
dt 



Elk 
T k 



(3) 



where r k is the relaxation time, assumed to be momen- 
tum dependent. Its momentum dependence will be de- 
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FIG. 1: (Color online) Dirac cones in graphene, in a situation 
of finite electronic density. The valence band is inert as long as 
the energy excitation processes are smaller than £ (the Fermi 
energy). 



tcrmined phenomenologically in such way that the de- 
pendence of the conductivity upon the electronic den- 
sity agrees with the experimental data. The Boltzmann 
equation is certainly not valid at the Dirac point, but 
since many experiments are done at finite carrier den- 
sity, controlled by an external gate voltage, we expect 
the Boltzmann equation to give reliable results when an 
appropriate form for r k is used. 

Let us compute the Boltzmann relaxation time, Tfe, for 
two different scattering potentials: (i) a delta function po- 
tential; (ii) a screened Coulomb potential. The relaxation 
time 77- is defined as (fe is the momentum) 



— = NtA 



k'dk' 



S(k, fc')(l ~ cos6*) , (4) 



where JVj is the total number of impurities in the sample, 
A is the area of the graphene sheet, and the transition 
rate S(k, fe') is given by 



S(k,k') 



2-7T 1 

-?-\H k ,, k \ 2 —6(k'-k), 
a vfH 



(5) 



where the vphk is the dispersion of Dirac fermions in 
graphene and H k ' k is defined as 



Hi 



fe' fe 



d 2 rr k -(r)U s (r)Mr), 



(6) 



with Us(i) the scattering potential and ipk(r) is the elec- 
tronic spinor wave function of a clean graphene sheet. If 
the potential is short ranged of the form Us = vo5(r), 
the Boltzmann relaxation time turns out to be 



4H 2 v F 1 



(7) 



where ni is the impurity concentration per unit area. On 
the other hand, if the potential is the screened Coulomb 
potential, given by Us(r) = eQe~ r / iD /(47re er) for 
charged impurities of charge Q, the relaxation time is 
given by 



1 

Tk 



Ak 2 L 2 D 



1 



h 2 k 



L y) 



(8) 



where «g 
obtains 



mQ 2 e 2 /(We 2 



In the limit Ljj —> oo one 



T k = 



(9) 



As we argue below, the phenomenology of Dirac fermions 
implies that the scattering in graphene must be of the 
form ^j. In what follows we explore the consequences of 
this type of relaxation time. 

At zero temperature the chemical potential £ (Fermi 
energy) is related to the density of charge carriers (see 
Fig. [T]) in the conduction band by n = ir^ 1 ((/vfH) 2 . 
In what follows we give the expressions for the transport 
coefficients in terms of £. 

The DC conductivity. Within the relaxation time ap- 
proximation the solution of the linearized Boltzmann 
equation when an electric field is applied to the sample 
is 



g k = 



df°(e k ) 
dtk 



eT k v k ■ E , 



and the electric current reads 



J = -r ^2 ev k9k ■ 



(10) 



(11) 



Since at low temperature the following relation 
— /°(efc)/<9efc — > <5(( — VF^k) holds, one can easily see 
that assuming Eq. © where k is measured relatively to 
the Dirac point, the electronic conductivity turns out to 
be 



2 a2 



e 2 tt(Hvf) 2 

i 9 



(12) 



where Uq is the strength of the scattering potential (with 
dimensions of energy). The electronic conductivity de- 
pends linearly on the particle density, in agreement with 
the experimental dataji^ and as was first noted by No- 
mura and MacDonald in Ref?i2,. The form used for Tk is 
thus imposed both by dimensional analysis and by phe- 
nomenology. We note that for short range scatterer o 17 i 18 
Tk oc fc" 1 , leading to an DC conductivity that is indepen- 
dent of the density and therefore in disagreement with the 
experiments. 

Having settled the need for Eq. ([9]) we now proceed to 
determine the density dependence of the other transport 
coefficients. We stress that the Coulomb potential is one 
possible mechanism of producing a scattering rate of the 



3 



form ([9]) but we do not exclude that other mechanisms 
may exist. 

The optical conductivity. Here we want to obtain the 
electronic density dependence of the optical conductiv- 
ity of a doped graphene plane. Since the Boltzmann 
approach does not include inter-band transitions, the ex- 
pressions obtained below are only valid as long as hw < ( , 
where the above mentioned transitions are blocked by the 
Pauli principle. 

Our aim is to obtain the response of the electronic 
system to an external electric field of the form 



E = E Q e 



i(q-r-ujt) 



(13) 



The Boltzmann equation has, for this problem, the form 

d/°(e fc ) „ 5fc . , 9g k f . 

S ev k ■ E = h v k ■ V r g h + — — . (14) 

de k r k at 

The solution of the linearized Boltzmann equation (Tl4]) 
is well known^ 2 -, reading 



de k 



with 



eT k v k ■ E Q 



ILUTk + ITkq ■ Vk 



(15) 



(16) 



The Fourier component, J(u>,q), of the current is given 
by 



J((J,q) = J d 2 kev k <S> q (uJ,k) 



df(e k ) 
dek 



, (17) 



leading in the long-wavelength limit to an optical con- 
ductivity of the form 



c 2 



■ iujhQ/u^ 



< 1 + (wfiC/ltg) 



(18) 



What should be stressed about Eq. (fT8|) is its density 
dependence. Since the electronic density is proportional 
to £ 2 it should be possible to collapse the optical conduc- 
tivity curves for different densities in a representation of 
o~xx(u)/n versus to^/n. This prediction can easily be con- 
firmed by reflectance measurements in graphene, for fre- 
quencies below £. Moreover, since the electronic density 
induced by a gate voltage V g is proportional to the gate 
voltage, by the relation n = ^P-Vg, where t is the thick- 
ness of the Silicon oxide substrate, one obtains a collapse 
of the optical conductivity for different gate voltages, by 
representing cr xx (u>)/V g versus u)y/V g . 

Thermal conductivity and thermopower. Our next goal 
is to obtain the density dependence of the thermal con- 
ductivity and of the thermopower. For a discussion 
within the Kubo-formalism, see Ref<». 

In the presence of a temperature gradient in the sam- 
ple, the linearized Boltzmann equation has the form 



df°(e k ) 
de k 



Vk' 



£fe - C 
T 



V r T + eE , 



bsv. 



9 ~, (19) 

Tfc 



where the measured electric field is given by E b sv , — 
E — V r £/e. In this situation we have, in addition to the 
electric current, a heat current (flux of heat per unit of 
area) given by 



4 x - 

U = -j2^Vk{ek - C)9k- 

k 



(20) 



Both the electric and the heat currents can be written 

„„22 



J = e 2 K -E obsv . + -K 1 -(-V r T) 

U = eK x ■ E obsv . + ^K 2 ■ (- V r T) , (21) 

where JQ, i = 0,1,2 are second order tensors. In this 
problem the tensors are diagonal, i.e. Ki = lfcj, and by 
a well established procedure 2 ^ one obtains 



fe 



2C 2 
h Uq 



(22) 
(23) 
(24) 



From the results J22J), ([23]) and (J24J) it is easy to derive 
both the thermal conductivity k and the thermopower 
Q. These are given by 



4tt 2 
3"^ 

27T 2 

3 h 



{k B Tf 



T { 3 h 



^T\:^BTf^~i^ik B Tf^ 



and 



eT3"C 



(k B Tf 



(25) 



(26) 



Again, what should be emphasized in these results is 
the dependence of both k and Q on the particle density, 
which is different from that of the usual two dimensional 
electron gas. Since it is experimentally feasible to control 
the carrier density in the graphene planei it is possible 
to check experimentally the dependence of the transport 
coefficients on the particle density. 

Hall effect. In what follows we prove that the scatter- 
ing rate ^ leads to a classical (low field) Hall coefficient 
R given by 



en 



(27) 



in agreement with the experiments. For a discussion 
within the Kubo-formalism, see Rcf. 24 . 

The linearized Boltzmann equation in the presence of 
a static magnetic field, perpendicular to the graphene 
plane H = (0, 0, B), is given by 



.^y)eE-v k = ^ + Uv k xH)-V hgh . (28) 
dek ) r k h 
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The determination of gj- becomes quite simple by realiz- 
ing that the right-hand-side of Eq. (|28p (after multipli- 
cation by Tfe) is the first order expansion of 

(29) 

Since Vf. = vpk/k and 77- is given by Eq. ([9]), Eq. 
(|29| can be rewritten as gk+(aev F /h)(kxH) , where a = 
vf1i 2 Uq 2 . We now define a new momentum variable k' 
given by 

h 

Solving Eq. (|30|) for k as function of fc' and making the 
left hand side of Eq. (f2"5)l equal to gw , one obtains the 
solution 

g(k) = (- d g^ J ev F ak-{E-aev F HxE/h)/(l+a 2 ) , 

(31) 

with a = Vpti 2 / (1 2 b Uq) and l 2 B = H/(eB). 

Knowing the solution to g(k) the electric current is 
easily computed giving 

J = -TT2 (E-azxE) , (32) 

where z is a unit vector along the z direction. The con- 
ductivity tensor <r is therefore given by 

a=-^( 1 °) (33) 
1 + a 2 V ~ a 1 / 

In the traditional Hall set up one has J y — leading to 
Jx = cr xx E x , (34) 



meaning that there is no magnetoresistance, and 



irv 2 h 2 1 
E y = -^—HJ X = —HJ X , (35) 



which produces a R Hall constant given by Eq. (f2"T|) . This 
result is in agreement with the experimental findings in 
the low field Hall effect. 

Conclusions. In this brief report we derived, using 
a semi-classical approach, the electronic density depen- 
dence of the DC conductivity, the optical conductivity, 
the thermal conductivity, the thermopower and the clas- 
sical (low field) Hall effect. Our proposed expressions are 
based on a phcnomcnological equation for the scattering 
rate of the Dirac electrons. Some of our findings have al- 
ready received experimental confirmation; the results for 
the thermal properties and for the optical conductivity 
can be tested experimentally. The Boltzmann approach 
can not explain the universal conductivity value occur- 
ring for n — > 0, that is, our approach breaks down at the 
neutrality point. 
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